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^-1 Abstract 

Some finite series of harmonic numbers involving certain reciprocals are eval- 
uated. Products of such reciprocals are expanded in a sum of the individual 
Lj reciprocals, leading to a computer program. A list of examples is provided. 
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1 Definitions and Basic Identities 

The generalized harmonic numbers used in this paper are: 

r- H^ = y -?- (i.i) 

V) fc=l 

in 

from which follows that Hq = 0. The traditional harmonic numbers are: 

^ m 

O H n = H™ (1.2) 



A well known identity is G3 |S] : 






Er jy * = ?( fl 2 + fl S ,) ) (°) 



and OH H]: 



and [5]: 



E rxi^-^fl^x-flSi) (1.4) 

fc=0 

n 

J2 \n n - k = El- i4 2) (1.5) 

fe=i 

n 

E fcTT^"- fe = H " +1 ~ ^ (L6) 



fc=0 



When a < b are two integers and {x k } and {y k } are two sequences of complex 
numbers, and {s k } the sequence of complex numbers defined by: 



k 

Sk = ^Xi (1.7) 



then there is the following summation by parts formula [1]: 

6-1 6-1 



^2 XkVk = Sb-iVb - ^2 s k (y k+1 - y k ) (1.8) 



k=a k—a 



2 Harmonic Number Identities with a Reciprocal 

Theorem 2.1. For nonnegative integer n and integer p > 0: 

" 1 1 

/ , jrp^ Hk = H n+p {H n+ i + H p -i) - -[(H n+ i + H p -if + H n+1 + Hp^] 

p-2 

n + k + 2 



k + p 
fe=o v 



E n + l . -^ 



fc=0 



Proof. Summation by parts (1.8) with x^ = l/(fc + p) and y k = i?fe yields: 



fc + 1 ^ fc + 1 (fc + l) s 



s=2 



Performing the summation over n and using (1.4) yields: 
" 1 1 

/I Trr^ Hk = H n+p H n+l - sC^n+l + #n+l) ~ -Hp + * 



fe =l fc+ P 2 



p-1 

s 



" 1 

E-(-Hn+«+l — #n+l — H s+ i + 1) 



(2-1) 



n -. n 1 

/ , k + v Hk = ( H n+P ~ H p )H n+ i - 2_^ -j^—{H k+p - Hp) (2.2) 



fc=i J fe=i 

Using: 

p 1 

# fe+p = #* + X! ^^ ( 2 - 3 ) 

s=l 

and for s > 1: 

1 = J_(^ L_) (24) 

(fe + s)(fc + l) s-l V fc + l fc + s ; V • ; 

yields: 

P 

^^^^(k^^^kh-kh^ (2 - 5) 



(2.6) 



s=l 



Using 



H, 



n+s + l — *l n+ 



H„ 



s 1 

_1 ^ n + k + 1 



(2.7) 



fc=i 



and changing the order of summation over s and k 

p— 1 s _, p— 1 -, p— l 



8=1 fc=l 



vly — — = y — — V- 

^ s ^ n + k + l *-^ n + k + 1 ^ s 

k—l s—k 

p-1 x 
f-^ n + k + 1 

P-2 | 

= ff p _i(iJ„ +p - ff„+l) - 2J — TT 



(2.8) 



fc=0 



n + fc + 2 



:H h 



and using ff s+1 = H s + l/(s + 1) and Q and l/(s(s + 1)) = 1/s - l/(s + 1) yields 
the theorem. □ 

Theorem 2.2. For nonnegative integer n and integer p > 0: 
n 1 1 



fc=0 



p-2 






(2.9) 



fc=0 



n + fc + 2 



-Hi: 



Proof. Summation by parts (1.8) with Xk = l/(?i + p — fc) and y^ = Hk yields: 



n -. n -. 

E^— #ra-& = 2, , 
A- -I- n z » 77 4- 77 



fc=0 



fc + p 



fe=l 



n + p — k 



Hi. 



n 1 

{H n+ p-l - H p _ 1 )H n+ i - 2_^ T , 7 {Hn+p-l 



(2.10) 



H 



n+p— £;— 1 J 



Using: 



and for s > 0: 



fc=i 



p-i 



H 



n+p— fc — 1 — ^n—k 



' J n -4- q 



1 



1 



n 4- 3 — k 



1 



1 



(n + s - k)(k + 1) n + s + 1 k + 1 n + s - k 



:) 



(2.11) 



(2.12) 



yields: 



fc + 1 " +p fel fc + 1 " fe ^n + s + rfc + 1 



1 

71 + S — fc 



) (2-13) 



Performing the summation over n and using (1.6) yields 

1 
k + p 



n ^ 
2 , u,„Hn-k = Hn+p-1 - H p -iH n+ i + H n+1 - H n+1 - H n 



fc=l 



(2.14) 



Using 



+ / J — ; — r{Hn+s-i + H n+ i — H a _\ — 1) 

'-^ n + s + 1 

s— 1 

p-1 J 

/ , ; —- 7 (Hn+1 - 1) = (-Hn+1 - 1)(-Hn+P _ H n+1 ) 

* — ' n + s + 1 



(2.15) 



and#„ +s _i = JT„ +s -l/(n + s) and l/((n + .s)(n + .s + 1)) = l/(n + s)-l/(n + S + l) 
and with ( 1.6 1: 



p— 2 n+p— 1 n -. 

/ . — ; 7~7;H n+s+ i = 2, i — ~~7H k — 2, -, — — 7 

*— ' n + s + 2 '-^ fe + 1 ^^ « + 1 



F fc 



s=0 



(2.16) 



I Til _ U\ 2 ) __ Zj2 i If\*) \ 

yj^n+p n n+p n n+l ~r n n+ll 



r(2) 



yields the theorem. □ 

Theorem 2.3. For nonnegative integer n and integer < p < n: 

" 1 1 1 

E ^^ - 2 [(jff - p+1 + Hp)2 + fr »-*+ 1 + ^ )] ~ ^"+ i(j ^ + n - p +i ) 



fc= P +i 



p-i 

E 

fc=0 



1 

n — p + k + 2 



:-fffc 



(2.17) 



Proof. Summation by parts (1.8) with Xk = l/(k — p) and y& = #& yields: 

n 1 n .. 

^ k~^o Hk = Hn + iHn -p - E k+i Hk - p 

F k=p+l 

n—p _. 

= H n+1 H n _ p - ]T fc+p + 1 #* 



k=p+l 



(2.18) 



fc=i 



The last sum is (2.1) with p replaced by p + 1 and n by n — p, which yields the 
theorem. □ 

Theorem 2.4. For nonnegative integer n and integer < p < n: 

n 

E u _ n Hn - k = H n-p - H n _ p (2-19) 



fc=p+l 



k — p 



Proof. 

n—p—l 



H 'k~ P Hn - k = ^ fc^Y jff "-p- fc - 1 ( 2 - 2 °) 



/c=p+l l k=0 



The last sum is (1.6) with n replaced by n — p — 1, which yields the theorem. □ 



3 Products of Reciprocals 

A finite product of these reciprocals with different p's can be written as a sum of the 
individual reciprocals. The formula for two reciprocals is, where p\ ^ p 2 : 



1 1 (/c+pi)-(fc+p 2 ) 



(fc+Pi)(fc+p 2 ) V\~V2 (k + Pi)(k+p 2 ) 
1 (^ l -) 



(3.1) 



Pi - P2 k+p 2 k+pi 
The formula for three reciprocals is, where Pi 7^ p 2 7^ Ps- 

1 t-V-^Vh ' 



(fc+Pl)(fc+p 2 )(fc+P3) Pl~P2 P2~PZ Pl~P3 k+p 3 

11 11, 



(3.2) 



P2-p-ik+p 2 pi - p 3 k+pi 

The recursion formula for to reciprocals in terms of the formula for to — 1 reciprocals 

is: 

m—X -. m— 1 -. 

2—1 2—1 

?72 1 rn— 1 1 .. 

TT j; = V^ 1 i 

M * + P* ^ l k + ip rn k+pi 

2—1 2=1 



m-1 

* — » n : — n W -4- Tt : 



j=l 



Pi - Pm k + Pi 



1 ^ 1 

+ 7— — 2^ a * — ~ — 

n T Pm . Pi Pm 

This recursion formula means that starting with rn — 1 and ai = 1, in each pass 
for certain to > 1 the a* for i = 1 ■ • • rn — 1 are divided by p m — pi, after which a m 
is minus the sum of the new o^ for i = 1 • • • m — 1. This way the recursion formula 
reduces to a double iteration, and it is also clear from this that for to > 1: 



J2 a * = ° ( 3 - 5 ) 



When the on have been computed, each individual reciprocal can be summed using 
the appropriate formula in the previous section, where the following substitutions are 
made: 

# <+> = H^\ + Y -. L- (3.6) 

n+p n+1 £^ ( n + Mm *■ > 

P-l , 

fc=0 v ; 

After these substitutions the coefficient of H„ +1 injhe formula for each individual 
reciprocal is identical, and therefore, by equation (3.5), the coefficient of H^ +1 in the 



resulting formula for m > 1 is zero, which means that for these products of these 
reciprocals only terms linear in harmonic numbers remain. 

4 Examples 

n 

Erri^'a^ 1 "^!) (4- 1 ) 

fc=0 

y — #fc = -{H 2 n+l - H&) + —H n+1 - — (4.2) 

<^k + 2 2 V n+1 n+1 ' n + 2 + n + 2 V ' 

k—0 

V X H - 1 (H* ff( 2 M i 2 " + 5 ff (" + l)(7n + 20) 

1- fc + 3^ ~ 2 (H " +1 " F » +l) + (n + 2)(n + 3) S ' l+1 " 4(n + 2)(n + 3) (4 ' 3) 



fc=0 



n 

J2iH k = l(H*+HW) (4.4) 

2- fc^2 Hfe - 2 ( ^" +1 + F ™+ l} " ( n -l)n(n + l) ffn+1 + 4n(n + 1) (4 ' 6) 



fe=3 



™ 1 1 

£ (fc + i)(fc + 2) Hfe = ^T2 " n^ Fn+1 (4 ' 8) 



fc=2 

•,2 



V^ X tr _ 9n^ + 5n-2 1 

2- (fc-l)(fe-2) fffe - 4n(n + l) " W^l n+1 (4 ' 10) 



! n + 1 

"-"n+l 



V" 1 TT = If7( 2 ) 

^fc(fc + l)(fc + 2) k 2 ™ +1 2(n + l)(n + 2)"" T± 2(n + 2) 
^ (k + l)k(k-l 



5n + 3 



fc=2 



(fc + l)fc(fc-l) 4(n+l) 2n(n+l) 



1 






>p i _ 2n 2 + 2n - 1 1 

^ k(k-l)(k-2) k ~ 4n(n + l) 2n(n - 1) " +1 



E 

fc=2 



1 



■ff fc — 



23n 2 + 57n + 28 



(k + 2)(k + l)k(k-l) 36(n + l)(n + 2) 



' U l W" (2) 



3n(n + l)(n + 2) 



" 1 

E (A . + i) fc ( fc _i) 



1 



r(2) 



(/,— l)A-(7,--l)(/c-2) Ffe e^^ 1 3(n - l)n(n + 1) 

2n 2 + l 



#, 



n+l 



12n(n + l) 



(4.11) 
(4.12) 
(4.13) 

(4.14) 



(4.15) 



n 

E 1 TT _ tt2 rr(2) 

. , i n n-k — J^n+1 n n+l 
fe=0 

n 

fe=0 

n 1 

E 1 tt _ tt2 rr(2) 

t, o-tln-k — n n+l n n+l 



(4.16) 

(4.17) 



fe=0 



fc + 3 



3n 2 + 7n-2 n + l 

2(n + 2)(n + 3) " +1 (n + 2)(n + 3) k ' ' 

n 1 

J] fc#n-fc = #n " ^ (4-19) 

fc=l 

f 1 ff _^2 ff (2) 2(2n + l) 2(3n 2 + 3n+l) 

^ 2 fc^i^ i - fe - ff - +i _ ff - +i _ ^ny F " +1 + n 2 (n+i) 2 (420) 



V^ X FT _ H-2 rr(2) _ 2 ( 3i 

2^k-2 n ~ k ~ n+1 n+l ln-1 



2 „ {2 ) 2(3n 2 - 1) 



fc=3 



(n — l)n(n + 1) 



H n +i 



n 1 

" 1 

£(*Ti)( 



- n ~ 1 ff 

fc — "T-Hn+1 ^ 

n+l 



(fc + l)(fc + 2)' 



#, 



n— fc 



n + 2 



2(6n 4 - 3n 2 + 1) 


(4.21) 


(n — l) 2 n 2 (n + l) 2 


n — 1 


(4.22) 


(n + l) 2 


-"n+l 


(4.23) 



V 1 H - n - 2 H (™-2)(2n + l) 

^ Uk^T) Hn - k ~ — Hn+1 ~ n^n + l) (424) 

y< 1 _ n - 3 (n-3)(3n 2 -l) 

j^(fc-l)(Jfe-2) n ~*~n-l "+ 1 (n-l) 2 n(n+l) l j 

V 1 rr _ " 2 + 3n-2 3n-l 

^fc(Jfe + l)(fc + 2) n ~ fc 4(n+l)(n + 2) n+1 4(n + l) 2 l J 

y^ 1 _ n 2 +n-4 zj in 3 + 3n 2 - 9n - 4 

A, (fe + l)fc(fc - 1) ™- fc_ 4n(n + l) n+1 4n 2 (n + l) 2 l j 

y 1 ff _ n 2 - rc - 4 (n 2 -2n-l)(5n 2 + 3n-4) 

2.^ k(h-^(k-o\ n ~ k w«.-ir n+1 4r«_i^„2r„ _i_-n l ^ 



fc=3 



k(k-l)(k-2) n ~ K 4n(n-l) + 4(n - l) 2 n 2 (n + 1) 

Av 1 _ n 3 + 3n 2 + 2n - 12 

^ (fc + 2)(fc + l)fc(fc-l) iJn_fe = 18n(n + l)(n + 2) Hn+1 

7n 3 + 18n 2 - 25n - 12 



(4.29) 



36n 2 (n + l) 2 

,A 1 n? - n - 12 

^ (fc + l)fc(fc-l)(fc-2) ff ™" fe _ 18(n-l)n(n + l) Hn+1 

9n 5 + 9n 4 - 41n 3 - 81n 2 + 32n + 24 
36(n- l) 2 n 2 (n + l) 2 



(4.30) 



5 Computer Program 

The Mathematical [9] program used to compute the expressions is given below. 

HarmNumPlus [p_,m_] :=HarmonicNumber [n+1 ,m] +Sum[l/(n+k) ~m,{k,2,p}] 
HarmNumMinus [p_,m_] : =HarmonicNumber [n+1 ,m] -Sum[l/(n-k+l) ~m,{k,0,p-l}] 
HarmSumPPos [p_,d_] :=Simplify[ 

HarmNumPlus [p , 1] (HarmonicNumber [n+1] +HarmonicNumber [p-1] ) 

-1/2 ( (HarmonicNumber [n+1] +HarmonicNumber [p-1] ) "2 

+HarmonicNumber [n+1 ,2] +HarmonicNumber [p-1 ,2] ) 

-Sum [l/(n+k+2) HarmonicNumber [k] ,{k,0,p-2}] 

-Sum [1/ (k+p) HarmonicNumber [k] , {k , , d-1}] ] 
HarmSumPNeg [p_ , d_] :=Simplify[ 

1/2 ( (HarmNumMinus [p , 1] +HarmonicNumber [p] ) ~2 

+HarmNumMinus [p , 2] +HarmonicNumber [p , 2] ) 

-HarmonicNumber [n+1] (HarmonicNumber [p] +l/(n-p+l) ) 

+Sum[l/(n-p+k+2) HarmonicNumber [k] ,{k,0,p-l}] 

-Sum [l/(k-p) HarmonicNumber [k] ,{k,p+l ,d-l}]] 
HarmSumP [p_ , d_] : =If [p<=0 , HarmSumPNeg [-p , d] , HarmSumPPos [p , d] ] 



HarmSumQPos [p_,d_] :=Simplif y [ 

HarmonicNumber [n+1] (HarmNumPlus [p , 1] -HarmonicNumber [p-1] ) 

-l/2(HarmonicNumber [n+1] ~2-HarmNumPlus [p, 1] ~1 

+HarmonicNumber [n+1 , 2] +HarmNumPlus [p , 2] ) 

-Sum [l/(n+k+2) HarmonicNumber [k] ,{k,0,p-2}] 

-Sum [1/ (k+p)HarmNumMinus [k+1 , 1] , {k, , d-1}] ] 
HarmSumQNeg [p_ , d_] :=Simplify[ 

HarmNumMinus [p+1 , 1] "2-HarmNumMinus [p+1 , 2] 

-Sum[l/(k-p)HarmNumMinus[k+l,l] ,{k,p+l ,d-l}] ] 
HarmSumQ [p_ , d_] : =If [p<=0 , HarmSumQNeg [-p , d] , HarmSumQPos [p , d] ] 
HarmTable [m_] :=Table [HarmonicNumber [n+1, i] ,{i,m}] 
HarmSumPQ[s_Integer,f_] :=Module [{d,u,t=HarmTable [2] } , 

d=If [s<=0,-s+l,0] ;u=Factor [Coefficient Arrays [f[s,d] ,t]] ; 

u[[l]]+Dot[u[[2]] ,t]+Dot[Dot[u[[3]] ,t] ,t]] 
HarmSumPQ [s_ , f _] : =Module [{f acs , d , u , f uns=0 , l=Length [s] , t=HarmTable [2] } , 

f acs=Table [0 , {l}] ; f acs [ [1] ] =1 ; Do [Do [f acs [ [j] ] /= (s [ [i] ] -s [ [j] ] ) ; 

facs[[i]]-=facs[[j]],{j,l,i-l}],{i,2,l}]; 

d=Min[s] ;d=If [d<=0,-d+l,0] ;Do [funs+=f acs [[i]]f [s [ [i] ] ,d] ,{i, 1,1}] ; 

u=Factor [Coef f icientArrays [funs , t] ] ; u [ [1] ] +Dot [u [ [2] ] , t] ] 
HarmonicSumP[s_] :=If [Length [s]==l , HarmSumPQ [s [ [1] ] ,HarmSumP] , 

HarmSumPQ [s , HarmSumP] ] 
HarmonicSumQ [s_] :=If [Length [s]==l .HarmSumPQ [s [ [1] ] , HarmSumQ] , 

HarmSumPQ [s , HarmSumQ] ] 

(* Compute some examples *) 
HarmonicSumP [3] //TraditionalForm 
HarmonicSumP [{2,1,0,-1}] //TraditionalForm 
HarmonicSumQ [{-2}] //TraditionalForm 
HarmonicSumQ [{0,-1 ,-2}] //TraditionalForm 
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